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General string-theoretic considerations suggest that four-dimensional large-A gauge theories 
should have dual descriptions in terms of two-dimensional conformal field theories. However, for 
non-supersymmetric confining theories such as pure Yang-Mills theory, a long-standing challenge 
has been to explicitly show that such dual descriptions actually exist. In this paper, we consider 
the large-A limit of four-dimensional pure Yang-Mills theory compactified on a three-sphere in the 
solvable limit where the sphere radius is small compared to the strong length scale, and demonstrate 
that the confined-phase spectrum of this gauge theory coincides with the spectrum of an irrational 
two-dimensional conformal field theory. 


Introduction. Gonfining gauge theories in the large- 
N limit are believed to have dual descriptions as weakly- 
coupled string theories [1] . Since string theories have 2D 
worldsheet conformal field theory (CFT) descriptions, it 
is expected that conhning 4D gauge theories may have 
alternative descriptions based on 2D CFTs. However, for 
non-supersymmetric quantum field theories (QFTs) such 
as Yang-Mills (YM) theory, no concrete relation between 
large-Af conhning theories and 2D CFTs has ever been 
found. 

In this paper we tackle this problem by studying the 
large-Af limit of 4D pure SU{N) YM theory, formulated 
at temperature T = and compactihed on a three- 
sphere of radius R. One can thus view the theory 
as living on x with Euclidean metric signature. 
The virtues of this setting are two-fold. First, thanks to 
asymptotic freedom, if we take AR ^ 1 where A is the 
YM strong scale, then the ’t Hooft coupling X = g^N 
becomes small — i.e. X{l/R) —0. As a result, the 
theory becomes solvable for any temperature j3 ~ . 

Second, it is known [2] that large-fV YM theory stays in 
the conhned phase when /3/i? > 1, even when A —)■ 0. In 
this context “confinement” means that the system has 
an unbroken center symmetry and that its free energy 
scales as iV°. As sketched in Fig. 1, it is plausible that 
the physics of YM theory is smooth as a function of AR. 
Thus, the AR 1 regime of the large-Af confined phase 
represents a particularly tractable 4D starting point in 
our search for a dual 2D description. 

Rather than attempt a string-theory construction of 
a 2D dual for large-A^ YM theory, we shall instead an¬ 
alyze the confined-phase spectrum of YM theory in the 
solvable AR 1 limit. Remarkably, we find that a sim¬ 
ple 2D CFT description emerges. Thus, in this limit, we 
conclude that the large-A^ confined-phase spectrum of 4D 
YM theory coincides with the spectrum of a 2D CFT. In 
the conclusions we briefly comment on possible relations 
between our result and string-theoretic expectations. 

Specifically, recall that the complete spectrum of a 
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FIG. 1. A conjectured phase diagram for large-A YM theory 
on S^ X S^. In the analytically tractable regime AR 1, the 
deconfinement transition occurs at ~ A, while for AR 3> 1, 
lattice studies have shown that it occurs at /3 ~ 1/A. This 
sketch illustrates the natural conjecture that these two limit¬ 
ing cases are smoothly connected. The results of this paper 
apply in the AR —^ 0 region indicated by the blue line. 

QFT is encoded in its grand-canonical thermal partition 
function. We take 4D YM theory to be minimally cou¬ 
pled to the metric, so that the Kaluza-Klein energies 
on the three-sphere are given by En = n/R in the A —>■ 0 
limit [2]. The partition function then takes the form 

OO OO 

Zym{P/R) = (1) 

n—0 n—0 

where q = and dn counts the number of states 

with energy En- Our main result is the demonstration 
that the grand-canonical partition function Zym of Yang- 
Mills theory coincides with a chiral partition function of 
a 2D CFT: 

Zym{t) = Z2d{t) ■ (2) 

In writing Eq. (2), we have analytically continued q to 
g27riT r € H, the complex upper half-plane. Thus 
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Im T = /?/ {2ttR) . Determining the physical interpreta¬ 
tion of Re r on the 4D gauge-theory side of Eq. (2) is an 
important problem for future work. 

The 4D partition function. We begin by briefly 
explaining the computation of Zym-, leaving a more 
leisurely exposition to Ref. [3]. To calculate the 4D par¬ 
tition function Zym('^)) we take the large-iV limit with 
A held fixed, which means taking the continuum limit 
after the large-limit. We work on x and as¬ 
sume that /3 and R are independent of N. Likewise, 
we do not consider states with energies > N because 
they lie beyond our UV cutoff. As is typical in studies 
of large-TV theories, we work with the U(N) version of 
YM theory rather than the SU{N) version [4]. When 
AR —>■ 0, the microscopic degrees of freedom of YM the¬ 
ory reduce to an infinite collection of color-adjoint-valued 
harmonic oscillators. These oscillators are counted by the 
massless-vector partition function, which can be written 
as 2 „(r) = — 2q^)/(l — q)^. The physical states 

are then determined by imposing the color Gauss law. 
In the A = 0 confined phase, the physical single-particle 
states can be identified with single-trace operators, and 
their energies are proportional to their scaling dimen¬ 
sions. The counting problem for these states, and also for 
the multi-particle states, has been solved [2, 5], and the 
resulting grand-canonical confined-phase partition func¬ 
tion is given by 


oo ^ 

^Ym(t) = ^ 


Hi 

n—1 

= l + 6q^ + I6q^ + 72q^ 


(1 - g”)3 


n—1 




— + q 


r,3n 


(3) 


Zym to be written as 


OO 

Zym = 

n—1 


(1 -I- ( 7")(1 — q'^z){l — q^z~^) 


(5) 


This observation is very important because the struc¬ 
ture of Eq. (5) matches the structure of the product 
representations of the Dedekind ry-function and gener¬ 
alized Jacobi d-functions. (In the related context of 
adjoint QCD, this was also noted in Ref. [ 8 ].) Specifi¬ 
cally, the Dedekind 77 -function has the product represen¬ 
tation ? 7 (t) = ~ while the generalized 

1 ?-funCti 0 n 7 T[ 3 ](r) = g("+“)V2g27r*n/3 ^ 

uct representation of the form 




n—1 


( 1 - 9 ") 


X (1 -h -k 


( 6 ) 


Under the S : t —>■ —1/t and T : t —>■ t + 1 genera¬ 
tors of the modular group SL{2,Z), we find 77 (—1/r) = 
y/—iTr]{T) and r]{T -f 1 ) = while 

5 : 7?[“](-1/t) = , 

T: dH(r + l)=e-“>[, + Ai/J(T) . (7) 


Given these definitions, the structure of Eq. (5) allows 
us to rewrite the 4D partition function Zym as a hnite 
product of Dedekind 77 -functions and Jacobi 7 ?-functions: 


As expected in any confining large-A^ theory, we find 
that the grow exponentially for large n. Thus, there 
are Hagedorn singularities in Zym(t). In Eq. (3), we 
find dn ^ and En ^ n for large n, with C = log(2 -|- 
^/3) « 1.317. This contrasts with the behaviors dn 
e'^ and En ^ y/n that would arise for a string theory 
with a flat target space. Of course, we are not in flat 
space: the spacetime curvature is ~ 1 /R, which is of the 
same scale as the effective string tension a' ~ 1 /i?^ that 
follows from our spectrum. The scaling properties of 
in Eq. (3) imply that the leading Hagedorn singularity of 
Zym(I3, fJ-i) is at Ph/R = C,^i = 0. Consequently, at 
77 ./ = 0 , there must be a phase transition to a deconfined 
phase at Ph- This is discussed in detail in Refs. [2, 6 ]. 

Modular symmetries. We now observe that the de¬ 
nominator in Eq. (3) can be factorized with roots that 
are inverses of each other: 

l-3q^-3q^^ + q^^ = {1 + q^){l - q^z)il - q^/z) (4) 

where z = 2 -|- -\/3- This pivotal algebraic observation 
was first made in Ref. [7] in the context of uncovering a 
subtle “temperature-reflection” symmetry for Zym- For 
our purposes, however, the key point is that this allows 




77(r)3 


/-y2e-"'’77(r)\ 

) 



( 8 ) 


where b = 7log(z)/27r « 0.21*, where 7 T 2 (t) = [^ 0 ^] 

and where the identity 2r]{2T)^ = 77 (T) 7 ? 2 ('r) has been 
used in passing from Eq. (5) to Eq. ( 8 ). The fact that 
b is imaginary is the reason the degeneracy factors dn in 
Eq. (3) grow as dn ~ e*"". The expression in Eq. ( 8 ) — 
and our interpretation of this expression in terms of spe¬ 
cific 2D CFTs, as discussed below — are the key results 
of our paper, with many striking consequences. 

Modularity versus dimensionality. The first in¬ 
teresting implication of Eq. ( 8 ) becomes apparent upon 
realizing that it is extremely unusual for the partition 
function of a 4D theory to be expressible as a finite prod¬ 
uct of modular functions, as in Eq. ( 8 ). (See, e.g., Ref. [9] 
for an early discussion along these lines.) The large-|r| 
behavior of a modular function is tied, through the S 
modular transformation, to its behavior near |t| = 0 . 
For example, the Dedekind 77 -function has the large-|r| 
expansion 77 (r) = — q + ...); the S transforma¬ 

tion then requires this function to behave at small |r| as 
77 (t) ^ exp[—77r/(12r)]/-\/~*T- Similar statements can be 
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made for the i?-functions. Thus, if a partition function 
can be written as a finite product of modular ry-functions 
and i?-functions, then it must have the leading behavior 


lim 

arg T—>-7r/2 


lim log -^modular ('^) 

\ t \^0 


aRiP (9) 


for a constant a. This amounts to the statement that 
logZ„ioduiar T as T —>■ oo. This is indeed the expected 
behavior for a 2D QFT. However, this is certainly not the 
expected behavior for a 4D QFT, for which we generically 
expect [10] 

logZ4D ^ as T —)• oo . (10) 


In this sense, 4D QFTs whose partition functions can be 
written in terms of modular functions behave as if they 
were 2D QFTs, since they follow Eq. (9) rather than 
Eq. (10). 

If we were to reverse the order of limits on the left side 
of Eq. (9) and take the jrj —>■ 0 limit with argr = tt/2, 
pure YM theory would follow the scaling in Eq. (10). 
Such a limit cannot be studied from Eq. (3) due to the 
Hagedorn singularities, and the physics is governed by 
the deconfined phase [2]. Eor Yang-Mills theory, Eq. (9) 
is thus valid only with the order of limits indicated. We 
note that in other theories such as adjoint QCD with pe¬ 
riodic boundary conditions for fermions, the Hagedorn 
singularities do not lie along argr = 7r/2 [8]; the two 
limits then commute and these theories exhibit 2D be¬ 
havior in the sense of Eq. (9) irrespective of the order of 
limits [3]. 

Vacuum energy. Another major consequence of 
Eq. (8) is that the modular properties of the rj- and t)- 
functions fix the vacuum energy Eym of our large-V YM 
theory to be zero. 

To see this, we first recall that if we write the q- 
series expansion of a modular function /(r) in the form 
/ = ^ thought of as the 

2D vacuum energy. Its value is fixed by the modular 
properties of / and tied to the values of a„. Were one 
to abitrarily shift A —> A -|- c, the modular properties 
of /(r) would be ruined because the ^-transformation 
would map q'^ = to thereby prevent¬ 

ing q’^/{ t) from transforming as a modular function. 

Next, we observe that the vacuum energy associated 
to the ry-function is 1/24, while 'd[i] has vacuum energy 
a^/2. Summing the vacuum energies of the individual 
modular functions in Eq. (8) we obtain a striking result: 

Eym = 0 . (II) 


Indeed, this is the only value consistent with the q- 
expansion for Zym given in Eq. (3), provided that Eym 
is calculated in a renormalization scheme which is consis¬ 
tent with the modular properties of Zym made evident 
in Eq. (8). This value, Eym = 0, coincides with the re¬ 
sult implied by T-reflection symmetry [7] , and also agrees 


with a direct evaluation of the sum over the confined- 
phase spectrum of finite-temperature large-A YM theory 
compactihed on as performed in Ref. [11]. 

CFT interpretation. The striking modular struc¬ 
ture of Eq. (8) suggests that the spectrum of our 4D YM 
theory coincides with that of a chiral (e.g., left- or right- 
moving) 2D GET. This motivates the central question 
we shall now explore for the rest of this paper: what is 
the 2D CFT which gives rise to Eq. (8) , and thus gives a 
2D description of 4D YM theory in the large-V limit? 

Unfortunately, we will not be able to give a complete 
answer to this question. The reason ultimately has to do 
with the fact that many distinct CETs can have coinci¬ 
dent spectra without being equivalent. They may differ, 
for example, in their correlation functions. In general, the 
most important aspects of a given 2D CFT are governed 
by its central charge (conformal anomaly) c and its spec¬ 
trum of operator conformal dimensions hi, i = 
where n is the number of so-called “primary” helds in 
the CFT. Along with the explicit traces over states, 
knowledge of c and the h^’s goes a long way in nailing 
down relevant aspects of the CFT such as its selection 
rules and correlation functions. But partition functions 
are only sensitive to the combinations = hi — c/24, 
rather than the values of c and hi individually. Conse¬ 
quently, without additional assumptions about the CFT 
in question (such as the assumption of unitarity, which 
would additionally tell us that min {hi} = 0), this repre¬ 
sents a fundamental limitation on our ability to specify 
a unique CFT. 

We will therefore answer a different but related ques¬ 
tion: do there exist any 2D CFTs to which our large-A 
YM theory is isospectraH Remarkably, we shall show 
that at least one such 2D CFT indeed exists. To see this, 
we first recall that a free c = 1 scalar CFT has a chi¬ 
ral spectrum whose trace is given by l/r](T), while the 
Z 2 orbifold of this CFT has a chiral sector whose trace 
is (277 (t)/i?2(t))^/^. Furthermore, the direct product of 
two copies of the c = —26 be ghost CFT has a chiral 
spectrum whose trace is given by ? 7 (t)^. Perhaps the 
most challenging to interpret is the remaining factor in 
Eq. (8), specifically 


V2e-^^^r]{T) 

llUir) 


However, this can be identified as the trace of the chiral 
{ e.g., left-moving) states in the vacuum sector of the c = 
2 bosonic f3'y ghost GET recently explored in Ref. [12]. 
This is a logarithmic CFT [13], and it has a U{1) con¬ 
served charge. Thus the vacuum-sector chiral partition 
function of the c = 2 f3'y CFT depends on the choice of a 
complex fugacity 2 ; = . To match with our expres¬ 

sions for YM theory, we set /r/3 = 2TTi b = — log(2 -|- ^/S). 

Putting this together, we therefore conclude that the 
expression in Eq. (8) can be viewed as the trace over 
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the chiral spectrum of a theory which is the direct prod¬ 
uct of five known CFTs, one of which is irrational. This 
then justifies the central claim of this paper in Eq. (2): 
there is indeed an irrational 2D CFT which is isospectral 
to the finite-temperature large-A^ 4D YM compactified 
on in the AR —^ 0 limit. Aside from explaining our 
observations concerning Eym and the small-| t| behav¬ 
ior of Zymi this fact has an intriguing further implica¬ 
tion. Two-dimensional CFTs have infinite-dimensional 
symmetries, which always include the Virasoro symme¬ 
try. Eq. (2) then suggests that large-YM theory has a 
hidden Virasoro symmetry. It would be very interesting 
to demonstrate this explicitly within YM theory. 

Primary operator spectrum. We now collect infor¬ 
mation concerning the spectrum of conformal dimensions 
corresponding to the primary fields of this tensor- 
product CFT. Our approach proceeds by determining 
the diagonal modular-invariant associated with the ex¬ 
pression in Eq. (8), and then computing the eigenvalues 
of the modular T operator to extract . 

We begin by defining the quantities 

^ _ -y2e-“^"r7(r)4 2r,{T) no^ 

where {m, n} are relatively prime integers (a relationship 
which we shall henceforth denote m T n), and P{k) = 
i(l -b (—1)^), fc S Z. Thus P{k) = 0,1 for odd or even fc, 
respectively. The set {Tjn,n} is a basis for a vector space 
over the field C with two key properties: it contains the 
“seed term” in Eq. (8), and it is the minimal set which is 
closed under the action of the 5'L(2,Z) modular group. 

The first property follows by noting that Tq^i{t) coin¬ 
cides with Eq. (8). The verification of the second prop¬ 
erty proceeds in two steps. First, it can be shown that, up 
to overall phases and extraneous factors of yf—ir, the S 
and T modular transformations map Tm,n to T-n,m and 
Tm,n+mi respectively. Second, we observe that if {m,n} 
are relatively prime, then {—n,m} and {m,n + m} are 
also relatively prime. Since all modular transformations 
can be generated by sequences of S and T, it then follows 
that the full modular “orbit” of our seed term Tgq is con¬ 
tained within the set of coprime integers {m,n}. Indeed, 
it is also possible to demonstrate [3] that the modular 
orbit actually covers all coprimes. 

As a result, the minimal “diagonal” modular-invariant 
generated from Eq. (8) is given by 

^diagonal = (Imrj^/^ ^ . (14) 

m_Ln 



FIG. 2. The numerical values of Eq. (14) with |m|, \n\ < 10, 
plotted within the unit-g disk. 


for an isolated set of points corresponding to the Hage- 
dorn singularities. The numerical values of ^diagonal on 
the interior of the unit-g disk are shown in Fig. 2. 

In order to extract the spectrum of effective conformal 
dimensions we now rewrite .^diagonal in a basis of 

eigenfunctions of the modular T : rr-l-1 operator. We 
do this because such eigenfunctions x(t) will have eigen¬ 
values exp[27ri/i|“®^] under T, allowing us to read off the 
values of (mod 1). Fortunately, constructing eigen¬ 
functions of the T-operator from linear combinations of 
the Tm „’s in Eq. (13) is relatively straightforward. Since 

r™.„(r + 1) = ^ (15) 

we see that any linear combination which includes Tm.,n 
must also include Tm,n+rm Tm,n+ 2 m-i and indeed all 
Tm,n+km where k G Z. Our T-invariant linear combi¬ 
nations can therefore be indexed by an arbitrary integer 
m and a second integer £ A m obeying 0 < £ < |m|. 
Hence T-eigenfunctions can be constructed analogously 
to Bloch eigenfunctions, by summing over all components 
Tm,e+km with k G Z with a Bloch phase a € [0,1) C K: 

XmAc = • (16) 

fcGZ 

It then follows that 


The appearance of the factor of (Imr)^/^ is standard 
when combining holomorphic and anti-holomorphic com¬ 
ponents, such as our Tm,n factors, each of which has mod¬ 
ular weight k = 3/2. It also ensures that .^diagonal is fully 
modular-invariant. Moreover, it can be verified numeri¬ 
cally that the infinite sum in Eq. (14) converges except 


Xm,l,a{T+l) = Xm.,e,a{T) 


(ett) 


where 


I (eff) 1 

o 


1 - P(m) 2 
-+TO 


\bf 


(17) 


(18) 
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One might wonder whether {xm,e,a} is the complete 
set of T-eigenfunctions. However, we have verified this by 
checking that summing over Xm,i,a reproduces Eq. (14): 

^diagonal = X! X! / (19) 

mGZ 0<^<|m| 

£_Lm 

This confirms that Eq. (18) is the desired set of effec¬ 
tive conformal dimensions (mod 1) of the primary opera¬ 
tors in our GET. The fact that these dimensions depend 
on q: — a continuous real variable — confirms that we are 
dealing with an irrational GET [14]. Our observations 
are consistent with the 2D logarithmic GET interpreta¬ 
tion discussed above, since it is known that logarithmic 
GETs typically have a continuously infinite number of 
primary operators [15]. 

Outlook. We have presented evidence that the con¬ 
fined phase of finite-temperature 4D non-supersymmetric 
large-iV pure Yang-Mills theory compactified on a three- 
sphere of radius R has a remarkable modular structure, 
as exposed by Eq. (8). This has many interesting con¬ 
sequences, such as the fact that this 4D gauge theory 
is isospectral to an irrational 2D GET in the AR —)■ 0 
limit, as summarized in Eq. (2). Moreover, as we shall 
demonstrate in a separate paper [3], modularity in the 
sense of Eq. (8) and isospectrality to 2D irrational GETs 
as in Eq. (2) turn out to be generic properties of large- 
confined-phase gauge theories with adjoint massless mat¬ 
ter in the A —>■ 0 limit. In Ref. [3] we shall also show that 
this structure is present in the large-iV limit of the A/” = 4 
superconformal index. 

It is not clear how easily the 4D-2D relation that we 
found, as summarized through Eq. (2), fits with standard 
string-theoretic expectations. Erom a string-theoretic 
perspective one might have expected that it would be the 
single-trace partition function — which can be thought of 
as representing the fluctuations of a single string — that 
would have a simple 2D GET description, assuming one 
is possible. It is less clear why the grand-canonical par¬ 
tition function which takes into account all multi¬ 

trace states and hence represents the fluctuations of an 
ensemble of many strings, should have a 2D GET de¬ 
scription. Erom this perspective, our result in Eq. (2) — 
and the analogous relations that we shall find in Ref. [3] 
for other, adjoint-matter gauge theories — are even more 
remarkable. 

Our results suggest a large number of interesting top¬ 
ics for future research. Obviously, it would be very inter¬ 
esting to understand whether Eq. (2) has an explanation 
within string theory, perhaps by making contact with the 
ideas in, e.g., Refs. [16]. It is also important to under¬ 
stand whether our large-fV 4D-2D spectral equivalence 
extends to correlation functions, and to explore how it 
is related to other known 4D-2D relations, such as those 
discussed in Refs. [17]. Note that unlike the 4D-2D rela¬ 
tions discussed in the context of supersymmetric indices 


(which by construction focus on a subset of states of the 
4D theory), our 4D-2D relation in Eq. (2) involves the 
full thermal partition function and hence concerns the 
entire spectrum of the 4D theory. 

Another interesting direction would be to develop an 
understanding of the modular structure of expressions 
like Eq. (8) directly from a 4D point of view, perhaps 
by making use of ideas from, e.(?.. Refs. [18]. Given re¬ 
cent progress in the understanding of the bulk duals of 
2D GETs (see, e.g., Ref. [19]), it is tempting to won¬ 
der whether our results may help to uncover the bulk 
dual of YM theory and of other non-supersymmetric 4D 
adjoint-matter theories. It would also be interesting to 
understand the extent to which the continuous spectrum 
of primary operators in the 2D theory suggested by our 
analysis has an interpretation in 4D YM theory. Einally, 
there remains the very important question of determin¬ 
ing how the modularity of Zym and our ensuing 4D-2D 
relation might evolve for A > 0. 
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